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Abstract 

We have carried out inclusive measurements of A(A) production in two-photon processes at TRIS- 
TAN. The mean y/s was 58 GeV and the integrated luminosity was 265 pb~^. Inclusive A(A) samples 
were obtained under such conditions as no-electron, anti-electron, and remnant-jet tags. The data 
were compared with theoretical calculations. The measured cross sections are two-times larger than 
the leading-order theoretical predictions, suggesting the necessity of next-to-leading-order Monte- 
Carlo generator. 
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1 Introduction 



Charm pair production in two-photon processes at ^/s=58 GcV have been reported in references j^, ||, 
|, |, |. AU of these papers have reported larger cross sections than that based on the lowest-order (LO) 
two-photon theory 0|, when the Drees-Grassie parametrization |^ in the resolved-photon processes 
was taken into account. In our previous reports ^ ^ we showed that the charm-pair cross sections 
approximately agreed with theory when we took into account the next-to-leading-order correction (NLO) 
to the hard interactions the lower charm-quark mass, and the intrinsic gluon Pt inside a photon |p^ , 
together with the Levy-Abramowicz-Charchula set-1 (LACl) parametrization |0 for the resolved-photon 
process. This explanation of excesses, however, must be confirmed by further analysis. 

In this article we present an analysis of inclusive A(A) production in the two-photon reaction using 
data taken with the TOPAZ detector at TRISTAN. As was demonstrated in our Kg analysis strange 
particles enhance the charm fraction in a sample. This analysis is, thus, expected to provide another 
handle for checking the above-mentioned theoretical arguments. Moreover, since A(A) is, being a strange 
baryon, favorably produced from gluon jets llj, |l3|, we expect to observe the NLO effects directly by 
measuring the inclusive A (A) production rates. 

In addition, the technique used to tag a remnant-jet using a forward calorimeter [Q, which was 
developed in the analysis can also be applied in this analysis. We have thus derived inclusive 
A (A) cross sections for both direct and resolved-photon processes separately. 



2 Event selection 

The data used in this analysis were obtained with the TOPAZ detector at the TRISTAN e+e~ collider, 
KEKflsl, p^ . The mean ^/s was 58 GeV and the integrated luminosity was 265 pb^^. A forward 
calorimeter (FCL), which covered 0.98 < | cos6\ < 0.998 (9 is the polar angle, i.e., the angle with respect 
to the electron beam), was installed during the course of the experiment. The FCL was made of bismuth 
germanate crystals (BGO), and was used to anti-tag the beam electrons (positrons) and to tag hadrons 
(remnant -jets) |0. The integrated luminosity of the data with the FCL detector was 241 pb^^. 

A description of our trigger system can be found in reference |p^ . The requirement for the charged- 
track trigger was two or more tracks with an opening angle greater than 45-90 degrees. The Pt threshold 
for charged particles was 0.3-0.7 GeV, which varied depending on the beam conditions. 

The event-selection criteria were as follows: there had to be three or more charged particles (Pt > 0.15 
GeV, I cos6\ <0.77), the invariant mass {Wvis) of visible particles (| cos 6*1 <0.77) had to be greater than 
2 GeV, the event- vertex position had to be consistent with the interaction point, and the visible energy 
had to be less than 25 GeV. In total, 280673 events were selected. 



3 Monte-Carlo simulation 

In order to estimate the acceptances and backgrounds in this analysis, we used the following Monte- 
Carlo simulation programs. Details concerning the event generation of direct as well as resolved-photon 
and vector meson dominance (VDM) processes can be found in references H, Here, we just 

note the following points. For cc generation, we used the current charm-quark mass of 1.3 GeV to 
calculate the cross sections for point-like processes and a constituent charm-quark mass of 1.6 GcV for 
the hadronization procedure, and made a next-to-leading order (NLO) correction by factorization, the 
details of which can be found in references H, |^ . Light-quark generation was carried out by using the 
lowest order (LO) formula with a P™™ cut of 2.5 GeV. We used the parton density functions by LACl 
[ ll| for the resolved-photon process, because the experimental data have favored this parametrization 
[|]|, |. JETSET6.3 111] was used for hadroniztion of parton system. Here, we must mention that we 
used the default values for the baryon production parameters, such as a P{qq)/ P{q) ratio of 0.10 and a 
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BMB ratio of 0.5. Generated events were processed through the standard TOPAZ detector simulation 
program ||2^ , in which hadron showers were simulated with an extended version of GHEISHA 7 . The 
simulation's handling of hadronic interactions with nuclei has been updated to fit various experimental 
cross sections. 

Using the above-mentioned Monte-Carlo simulations, the trigger efficiency for the sum of the direct 
and resolved-photon processes was estimated to be 79%, 97% of which represented charged trigger events. 
The event-selection efficiency after the trigger was obtained to be 80%. 

To generate single-photon-exchange hadronic events we used JETSET6.3||l^ with the parameter val- 
ues given in reference 20 1. 



4 Tagging conditions 

The tagging conditions were as follows. For anti-electron tagging there had to be no energy deposit of 
more than QAEi, in |cos0| < 0.998 (anti-electron tag or anti-tag), where Eb is the beam energy. The 
rejected sample was called the electron tag sample. These selected events were from collisions of almost- 
real photons, for which the equivalent photon approximation was expected to be accurate at the 1% level. 
When the energy cut was lowered, the number of mistakenly rejected beam remnant hadron (remnant- 
jets) events became significant, as predicted by Monte-Carlo simulations. This implies the possibility 
of tagging the resolved-photon process by requiring, for instance, 500 MeV < EpcL < 0.25£'f,, where 
EpCL is the energy deposit in the FCL (remnant-jet tag or rem-tag). The yield of the remnant-jet tag 
events agreed with our Monte-Carlo simulation within 5% [^. An event selection without these two tags 
is hereafter called, a "no-electron tag" or a "no-tag" . Also, data without using the FCL detector were 
included in this sample. The fractions of anti-electron and remnant -jet tag events to no-tag events were 
obtained to be 97.6 and 47%, respectively. 

In the Monte-Carlo simulations we used the equivalent photon approximation with the photon fiux 
expression 

/7/e(^7) = ^(1 + (1-X,)2)ln-^"-- ^errA-X, 



where 



27rT ^ ^ " ^ 
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We set the P^^x limit at the smaller of P^^ + and the anti-tag limit [2-E^(l — x^){\ — cosBmax), 
where = 0.4, 9max = 3.2 degrees], where PT,q and ruq are the transverse momentum and the mass of 
a quark, respectively p^ . 

The tagging efficiency of the remnant jet for the resolved-photon process was estimated to be 76 ±3%, 
taking into account an estimate of the FCL noise (described later) . We tried to generate remnant partons 
using two techniques: one was along the beam direction; the other used a Gaussian distribution of the 
PT-width (0.44 GeV) with respect to the beam axis. These two methods differed in acceptance by only 
1%. Further, the tagging efficiency of the remnant-jet tag for the direct process was estimated to be 
20 ± 2% with the estimated level of FCL noise. 



5 A (A) inclusive analysis 

The charged-track selection criteria for the A(A) analysis were as follows: for each track Pt had to be 
greater than 0.15 GeV, | cos6\ had to be less than 0.77, and the closest approach to the interaction point 
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Figure 1: Invariant-mass spectra of p7r 's and their charge conjugations. The sohd histogram is for the 
no-electron tag, the dashed one is for the anti-electron tag, and the dotted one is for the remnant -jet tag. 



in the XY-plane (perpendicular to the beam axis) had to be greater than 0.5 cm. Using these selected 
tracks, we looked for opposite-sign pairs with an opening angle of less than 90 degrees, and carried out 
secondary vertex reconstructions three-dimensionally. Here, the dE/dx of one of two tracks had to be 
consistent with the proton assumption (Xp(p) < 10)) ^ind that of the other had to be consistent with the 
pion assumption (x^± < 10). We demanded that these pairs be consistent with the assumption that 
they came from event vertices with flight lengths greater than 3 cm. Finally, we rejected pairs which 
had vertices near (±3cm) to the inner pressure vessel {Rxy ~30cm) or the field cage {Rxy ^33cm) of 
the TPC [l^ in order to reduce the background from nuclear interactions; by generating events without 
any A(A)'s, we found that there were no significant contribution from such fake pairs after this cut. We 
also required that | cosd\ of the pair had to be less than 0.77. The invariant-mass distributions of these 
candidate pTr~ pairs and their charge-conjugated states (CC) are plotted in Figure || for the three tagging 
conditions, respectively. These invariant-mass spectra were fitted with the sum of a fifth-order polynomial 
and a Gaussian distribution, and the peak entries were obtained for no-, anti-, and remnant-jet tags to 
be 255 ± 20, 215 ± 19, and 123 ± 13 A(A)'s, respectively, in the Pt range between 0.75 - 2.75 GeV. The 
peak position and the width were consistent with the detector simulation. Among them, 2-AA, 0-AA, 
and 0-AA pairs were found. In order to derive the differential Pt cross sections we divided Pt into four 
bins, as shown in Table |^. 



6 Background subtractions 

The single-photon-exchange process produced a large background, especially for high-Py A(A)'s. This 
background could have been reduced by applying a cut on the total visible energy. However, we avoided 
this in order to keep the acceptance for high- Pt A(A)'s. Instead, the contamination from the single- 
photon-exchange process was estimated and subtracted using a Monte-Carlo simulation on a bin-by-bin 
basis. The background fractions for the no-tag sample were 3± 1, 5± 1, 12±3, and 16±3%, respectively, 
for the Pt bins shown in Table |l|; they were strongly Pt dependent. We also estimated these fractions for 
anti- and remnant-jet tags. They were consistent with the above-mentioned values within the statistical 
errors. 
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tag cond. 


no-tag 


anti-tag 


rem-tag 


rem-tag 


anti-rem-tag 


VDM subt. 






yes 


no 


no 


Pt range 






cross sections 




(GeV) 






(pb/GeV) 






0.75-1 


108± 30 


102± 30 


29± 11 


85± 28 


73± 28 


1-1.25 


46± 12 


42± 12 


17± 6 


35± 11 


25± 12 


1.25-1.5 


13. 2± 4.9 


12. 1± 4.7 


5.2± 3.1 


7.4± 4.1 


6.9± 5.2 


1.5-2.75 


5.1± 1.4 


3.5± 1.2 


2.3± 1.0 


2.4± 1.0 


1.2± 1.4 



Table 1: Differential cross section of A(A) versus Pt (GeV) [da/dPr (pb/GeV)], for | cos 6*1 < 0.77. Six 
cases are listed: no-tag , anti-electron tag, remnant-jet tag with VDM subtraction, anti-remnant-jet tag 
without VDM subtraction, and remnant-jet tag without VDM subtraction, which are described in the 
text. 

The background from beam-gas interactions was estimated using the off-vertex events in the beam 
direction. The beam-gas contribution for a no-tag sample was 22 ±3% on the average, and was subtracted 
from the data. This was mostly due to a vacuum leak in the beam pipe during some period. It was larger 
than in the case of the analysis For the anti- and remnant-jet tag samples, the above values 
became 26 ± 4 and 28 ± 5%, respectively. Without the anti-tag condition, the beam-gas contribution 
became slightly larger, suggesting that the electron-tagged events were cleaner. 

FCL noise hits were studied by analyzing random-trigger and Bhabha events. The probability of noise 
hits with EpcL >0.5 GeV was estimated to be 13.7%, while for hits with EpcL > OAEi, it reduced to 
be 0.1%. The FCL noise was also related to a vacuum leak. In the Monte-Carlo simulations we added 
noise hits randomly in accordance with the observed noise-hit probability in order to reliably estimate 
the tagging efficiencies. 

7 Systematic errors 

The systematic errors for the cross sections were estimated, bin by bin, as follows. For the trigger, we 
added some extra noise hits in the tracking chambers in the simulations. For the event selection and 
the A(A) reconstruction we changed the cut values and evaluated the systematic errors as cross-section 
differences. We also changed the pulse-height threshold in the TPC simulation in order to evaluate the 
effects on its tracking efficiency. We added the obtained systematic errors quadratically on a bin-by-bin 
basis. 

In order to investigate the effect of nuclear interactions in the material in front of the TPC, we 
compared the yields of A and A in the experimental data as well as in the Monte-Carlo data. The ratio 
[iV(A)/iV(A)] was 1.4 ± 0.3 in the experiment, while the Monte-carlo simulation predicted this value to 
be 1.43 ± 0.09, being consistent with the data. The deviation from 1.0 was considered to be due to 
an inelastic scattering of A with the material in front of the TPC. This occurred mostly in the low- Pt 
regions. The effect of nuclear interactions was corrected using the Monte-Carlo simulations. The total 
systematic errors were 16~21%, depending on Pt, of which the cut dependence in the event selection 
was the dominant source. These systematic errors were quadratically added to the statistical errors. 

We also checked the acceptance ambiguity due to the parametrization dependence of the resolved- 
photon processes by comparing the LACl jl^ and Drees-Grassie [DC] |Q parametrizations. The ac- 
ceptance difference was estimated to be 3.3%, which is small compared to the systematic errors given 
above. 
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8 Results 



The Pt difFcrcntial cross sections were obtained from the number of reconstructed A(A)'s in each bin, and 
its corresponding efficiency was estimated using the previously described Monte-Carlo simulations. They 
are listed in Table and plotted in Figures § (a) - (d) for the three tagging conditions and two subtraction 
schemes discussed afterward, respectively. Figure y (a) is for anti-electron tag events. In the remnant- 
jet tag events, the Monte-Carlo simulation predicted a significant contamination from VDM events. The 
tagging efficiency for the VDM process was estimated to be 65 ±6%, slightly less than that of the resolved- 
photon process. In addition, there was a large ambiguity in the cross section of the VDM process. We 
therefore calculated the cross sections using two subtraction schemes. Figure § (b) was obtained by 
subtracting the VDM contribution predicted by the Monte-Carlo simulation (VDM subtraction) for the 
remnant jet. Here, the "anti-remnant-jet tag" cross section was obtained by subtracting the remnant-jet 
tag cross section from that of the anti-electron tag. We do not show the cross sections for the anti- 
remnant-jet tag with VDM subtraction because of low statistics. Figures I (c) and (d) were obtained 
without VDM subtraction. In Figures || (b) and (c), the contribution of the direct process was subtracted. 
The histograms in Figures ^ (a) - (d) are the Monte-Carlo predictions: the cross-hatched, singly-hatched, 
and open areas are predictions for the direct, resolved-photon (LACl), and VDM processes, respectively. 



9 Discussions 

The fraction of charm events was studied using the above-mentioned Monte-Carlo simulations. We found 
that 70% of the events with Pt(A, A) > 1.5 GeV were of charm origin. On the other hand, only 35% of 
the events with charged tracks of Pt > 1.5 GeV were from cc pairs. Also, the Monte-Carlo simulations 
predicted that 55% of these high-P^ charm events originated from the direct process. The contribution 
of the resolved-photon process is higher than that in the Kg analysis ||^ . This is considered to be due to 
the existence of gluon jets in the resolved-photon process, even in the LO calculation. In this study we 
derived six types of cross sections using different tagging conditions and subtraction schemes. We can 
therefore separately compare each cross section with the theoretical prediction for each process. 

Firstly, about 30% of the high-Py (Pt > 1-5 GeV) events can be explained as electron-tagged events 
(see Table |]). Secondly, the anti-tag cross section [Figure || (a)] is two-times larger than the theoretical 
model prediction, though the spectrum shape seems to be consistent with it. Here, we must mention that 
this theoretical model reasonably explained the Kg inclusive cross section in the previous study . 

The cross sections with the remnant-jet and anti-remnant-jet tags also show discrepancies having a 
factor of ^2 compared to predictions of the direct, resolved, and VDM processes [Figures ^ (b)-(d)]. 

The VDM model is considered to be ambiguous, especially in predicting the total cross sections. 
However, the A(A)-spectrum of this process is softer than the others and also than the experimental 
data. We therefore cannot fit it to the experimental data by changing the normalization. Moreover, 
increasing the normalization factor makes the Kg cross-section discrepancy in the low-P^ region larger, 
and inconsistent with the experimental observation [^. 

In order to check whether the parton-density functions have anything to do with the discrepancy, we 
compared our remnant -jet-tag data [Figure (c)] with the predictions from six sets of parametrizations 



by Hagiwara, Tanaka, Watanabe, and Izubuchi [WHITl-6] |23 . A systematic analysis on the gluon 
distributions can be carried out using these parametrizations. We selected those parametrizations which 
showed higher charm cross sections, i.e., WHIT-1 and -4. The results are shown in Figures I (a) and (b), 
where the histograms are the predictions by the WHIT-1 and -4 parametrizations with P™™'s of 2.0 and 
2.5 GeV. There are some possible combinations (for example, WHIT-4 with P^*" = 2.0 GeV seems to 
fit the data better). Notice that these operations do not solve any discrepancies which appeared in the 
anti-remnant-jet-tag data. 
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1 1.5 2 2.5 1 1.5 2 2.5 

Pt(A,A) GeV 

Figure 2: Differential cross section of A(A) versus Pt (GeV) [da/dPr (pb/GeV)], for | cos 6*1 < 0.77. 
Five cases are plotted: (a) anti-electron tag, (b) remnant-jet tag with the VDM and direct process 
subtraction, (c) remnant -jet tag without the VDM subtraction, and (d) anti-remnant -jet tag without the 
VDM subtraction, as described in the text. The corresponding processes arc presumably; (b) resolved- 
photon, (c) resolved-photon and VDM, and (d) VDM and direct processes. The histograms are the 
theoretical predictions which are described in the text. The open area is for the VDM, the singly-hatched 
one is for the resolved-photon process, and the cross-hatched one is for the direct process. 
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Figure 3: Differential cross section of A(A) versus Pt (GeV) [P^da/dPT (pb-GeV^)], for | cos 6*1 < 0.77 
for the remnant-jet tag without the VDM subtraction. The hatched areas are the predictions by the 
VDM Monte-Carlo simulation. The histograms are predictions by the WHIT-1 and -4 parametrizations. 
The solid one is WHIT 1 and the dotted one WHIT-4. Two values of P^™'s were used, i.e., (a) 2.0 and 
(b) 2.5 GeV. 
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tag cond. 


Experiment 


Theory (LO) 


Exp. /Theory 


subprocess 


antitag 


43.3±8.3 


19.1 


2.26±0.43 


VDM+resolved+direct 


rem-tag (-VDM) 


15.6±3.5 


6.0 


2.60±0.58 


resolved 


rem-tag 


34.8±7.8 


17.3 


2.01±0.45 


VDM+resolved 


anti-rem 


27.7±7.9 


13.1 


2.11±0.60 


VDM+direct 



Table 2: Total cross section (pb) of A(A) in the | cos6'| < 0.77 and 0.75 < Pt < 2.75 GeV range. The 
notation (-VDM) means the VDM subtraction which was described in the text. Here, we use the LO 
theories in order to show the discrepancy with the experimental data. 



In order to reduce these discrepancies we may need an NLO correction to the light-quark events, 
which is absent from our present Monte-Carlo generator. 

Our Monte-Carlo generation of 77 events was based on the LO matrix element calculations in the 
hard processes. Although the total cross section of the direct process was increased by a factor of 
1.31 & |[ ||, and the Pt{Pt of parton)-dependent factor was used in the resolved-photon processes 

Hi^i our Monte-Carlo generator has no explicit hard gluon emissions, which should exist in a truly 
NLO generator. We consider that the discrepancies between the data and the theory can be attributed 
to these effects, based on the baryon abundance in the experimentally observed gluon jets ||l^, |l^. The 
baryon excess in the gluon jets has not yet been established in a sufficiently quantitative way. When it 
is, our data will be direct evidence of the NLO effect. 

We show in Table || the total cross sections of A(A) in the | cos 6*1 < 0.77 and 0.75 < Pt < 2.75 GeV 
range. In order to extract the NLO effect, we calculated the LO cross sections, which are also shown in 
Table |[ Here, we removed the NLO factorization from our Monte-Carlo generator. The ratio between 
the theoretical predictions and the experimental data are typically ~2. 



10 Conclusion 

We carried out an inclusive measurement of A(A) productions in two-photon processes at TRISTAN. 
The mean y/s was 58 GeV and the integrated luminosity was 265 pb^-'^. Inclusive A(A) cross sections 
in two-photon processes were obtained under such conditions as no-electron, anti-electron, and remnant- 
jet tags. In particular, using the remnant-jet tagging we could unambiguously extract the contribution 
from the resolved-photon process. Comparisons with theoretical predictions were carried out. The cross 
sections were two-times larger than the theoretical predictions without hard gluon emissions, suggesting 
the necessity of next-to-leading-order Monte-Carlo generator. 
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